RT-PCR genotyping assays to identify SARS-CoV-2 variants in England in 2021: a design and retrospective evaluation study.

BACKGROUND: Whole-genome sequencing (WGS) is the gold standard diagnostic tool to identify and genetically characterise emerging pathogen mutations (variants), but cost, capacity, and timeliness limit its use when large populations need rapidly assessing. We assessed the potential of genotyping assays to provide accurate and timely variant information at scale by retrospectively examining surveillance for SARS-CoV-2 variants in England between March and September, 2021, when genotyping assays were used widely for variant detection. METHODS: We chose a panel of four RT-PCR genotyping assays to detect circulating variants of SARS-COV-2 in England and developed a decision algorithm to assign a probable SARS-CoV-2 variant to samples using the assay results. We extracted surveillance data from the UK Health Security Agency databases for 115 934 SARS-CoV-2-positive samples (March 1-Sept 6, 2021) when variant information was available from both genotyping and WGS. By comparing the genotyping and WGS variant result, we calculated accuracy metrics (ie, sensitivity, specificity, and positive predictive value [PPV]) and the time difference between the sample collection date and the availability of variant information. We assessed the number of samples with a variant assigned from genotyping or WGS, or both, over time. FINDINGS: Genotyping and an initial decision algorithm (April 10-May 11, 2021 data) were accurate for key variant assignment: sensitivities and PPVs were 0·99 (95% CI 0·99-0·99) for the alpha, 1·00 (1·00-1·00) for the beta, and 0·91 (0·80-1·00) for the gamma variants; specificities were 0·97 (0·96-0·98), 1·00 (1·00-1·00), and 1·00 (1·00-1·00), respectively. A subsequent decision algorithm over a longer time period (May 27-Sept 6, 2021 data) remained accurate for key variant assignment: sensitivities were 0·91 (95% CI 0·74-1·00) for the beta, 0·98 (0·98-0·99) for the delta, and 0·93 (0·81-1·00) for the gamma variants; specificities were 1·00 (1·00-1·00), 0·96 (0·96-0·97), and 1·00 (1·00-1·00), respectively; and PPVs were 0·83 (0·62-1·00), 1·00 (1·00-1·00), and 0·78 (0·59-0·97), respectively. Genotyping produced variant information a median of 3 days (IQR 2-4) after the sample collection date, which was faster than with WGS (9 days [8-11]). The flexibility of genotyping enabled a nine-times increase in the quantity of samples tested for variants by this method (from 5000 to 45 000). INTERPRETATION: RT-PCR genotyping assays are suitable for high-throughput variant surveillance and could complement WGS, enabling larger scale testing for known variants and timelier results, with important implications for effective public health responses and disease control globally, especially in settings with low WGS capacity. However, the choice of panels of RT-PCR assays is highly dependent on database information on circulating variants generated by WGS, which could limit the use of genotyping assays when new variants are emerging and spreading rapidly. FUNDING: UK Health Security Agency and National Institute for Health Research Health Protection Research Unit in Emergency Preparedness and Response.

Preliminary investigation of a significant national Cryptosporidium exceedance in the United Kingdom, August 2023 and ongoing.

Routine laboratory surveillance has identified an unprecedented and ongoing exceedance of Cryptosporidium spp. across the United Kingdom, notably driven by C. hominis transmission, since 14 August 2023. Information from 477 reported cases in England and Wales, followed up with a standardised exposure questionnaire as of 25 September 2023, identified foreign travel in 250 (54%) of 463 respondents and swimming in 234 (66%) of 353 cases. A significant, common exposure has not yet been identified in first analyses.

Exploring the Spatial Relative Risk of COVID-19 in Berlin-Neukölln.

Identifying areas with high and low infection rates can provide important etiological clues. Usually, areas with high and low infection rates are identified by aggregating epidemiological data into geographical units, such as administrative areas. This assumes that the distribution of population numbers, infection rates, and resulting risks is constant across space. This assumption is, however, often false and is commonly known as the modifiable area unit problem. This article develops a spatial relative risk surface by using kernel density estimation to identify statistically significant areas of high risk by comparing the spatial distribution of address-level COVID-19 cases and the underlying population at risk in Berlin-Neukölln. Our findings show that there are varying areas of statistically significant high and low risk that straddle administrative boundaries. The findings of this exploratory analysis further highlight topics such as, e.g., Why were mostly affluent areas affected during the first wave? What lessons can be learned from areas with low infection rates? How important are built structures as drivers of COVID-19? How large is the effect of the socio-economic situation on COVID-19 infections? We conclude that it is of great importance to provide access to and analyse fine-resolution data to be able to understand the spread of the disease and address tailored health measures in urban settings.

Animal and environmental risk factors for sporadic Shiga toxin-producing Escherichia coli (STEC) infection in England: a case control study for O157, O26 and other STEC serotypes.

Most Shiga toxin-producing E. coli (STEC) infections are sporadic. Routine enhanced surveillance questionnaires of confirmed STEC cases in England contained promising data to conduct a case-control study to identify non-food exposures linked to the risk of becoming infected with different STEC serotypes, including O157, O26 and all others; this study pulled eligible cases from the recorded enhanced surveillance data. Controls were recruited from the general population and answered a comparable postal questionnaire. Logistic regression was performed to identify risk factors associated with STEC infection for O157, O26 and other serotype cases. In adjusted models, travel outside of the U.K. and childcare occupations raised the risk of infection for all serotypes. Day trips within the UK, exposure to dogs and contact with soil were linked to lower infection risk. Resident region within England was often linked to decreased risk. Summer season was linked to O157 and O26, but not other STEC. Swimming in the sea was linked to increased risk of infection by O157, but not other types of STEC. Correlations between exposures and infection were similar when the analysis was repeated excluding participants with a history of foreign travel. As the first case-control study in England to include sporadic non-O157 STEC, the varying risk factors between O157 and non-O157 cases suggest there are potentially unique reservoirs for different serotypes.

Enhancing epidemiological surveillance of the emergence of the SARS-CoV-2 Omicron variant using spike gene target failure data, England, 15 November to 31 December 2021.

When SARS-CoV-2 Omicron emerged in 2021, S gene target failure enabled differentiation between Omicron and the dominant Delta variant. In England, where S gene target surveillance (SGTS) was already established, this led to rapid identification (within ca 3 days of sample collection) of possible Omicron cases, alongside real-time surveillance and modelling of Omicron growth. SGTS was key to public health action (including case identification and incident management), and we share applied insights on how and when to use SGTS.

Covid-19 Vaccine Effectiveness against the Omicron (B.1.1.529) Variant.

BACKGROUND: A rapid increase in coronavirus disease 2019 (Covid-19) cases due to the omicron (B.1.1.529) variant of severe acute respiratory syndrome coronavirus 2 in highly vaccinated populations has aroused concerns about the effectiveness of current vaccines. METHODS: We used a test-negative case-control design to estimate vaccine effectiveness against symptomatic disease caused by the omicron and delta (B.1.617.2) variants in England. Vaccine effectiveness was calculated after primary immunization with two doses of BNT162b2 (Pfizer-BioNTech), ChAdOx1 nCoV-19 (AstraZeneca), or mRNA-1273 (Moderna) vaccine and after a booster dose of BNT162b2, ChAdOx1 nCoV-19, or mRNA-1273. RESULTS: Between November 27, 2021, and January 12, 2022, a total of 886,774 eligible persons infected with the omicron variant, 204,154 eligible persons infected with the delta variant, and 1,572,621 eligible test-negative controls were identified. At all time points investigated and for all combinations of primary course and booster vaccines, vaccine effectiveness against symptomatic disease was higher for the delta variant than for the omicron variant. No effect against the omicron variant was noted from 20 weeks after two ChAdOx1 nCoV-19 doses, whereas vaccine effectiveness after two BNT162b2 doses was 65.5% (95% confidence interval [CI], 63.9 to 67.0) at 2 to 4 weeks, dropping to 8.8% (95% CI, 7.0 to 10.5) at 25 or more weeks. Among ChAdOx1 nCoV-19 primary course recipients, vaccine effectiveness increased to 62.4% (95% CI, 61.8 to 63.0) at 2 to 4 weeks after a BNT162b2 booster before decreasing to 39.6% (95% CI, 38.0 to 41.1) at 10 or more weeks. Among BNT162b2 primary course recipients, vaccine effectiveness increased to 67.2% (95% CI, 66.5 to 67.8) at 2 to 4 weeks after a BNT162b2 booster before declining to 45.7% (95% CI, 44.7 to 46.7) at 10 or more weeks. Vaccine effectiveness after a ChAdOx1 nCoV-19 primary course increased to 70.1% (95% CI, 69.5 to 70.7) at 2 to 4 weeks after an mRNA-1273 booster and decreased to 60.9% (95% CI, 59.7 to 62.1) at 5 to 9 weeks. After a BNT162b2 primary course, the mRNA-1273 booster increased vaccine effectiveness to 73.9% (95% CI, 73.1 to 74.6) at 2 to 4 weeks; vaccine effectiveness fell to 64.4% (95% CI, 62.6 to 66.1) at 5 to 9 weeks. CONCLUSIONS: Primary immunization with two doses of ChAdOx1 nCoV-19 or BNT162b2 vaccine provided limited protection against symptomatic disease caused by the omicron variant. A BNT162b2 or mRNA-1273 booster after either the ChAdOx1 nCoV-19 or BNT162b2 primary course substantially increased protection, but that protection waned over time. (Funded by the U.K. Health Security Agency.).

Evidence of on-going transmission of Shiga toxin-producing Escherichia coli O157:H7 following a foodborne outbreak.

In August 2019, public health surveillance systems in Scotland and England identified seven, geographically dispersed cases infected with the same strain (defined as isolates that fell within the same five single nucleotide polymorphism single linage cluster) of Shiga toxin-producing Escherichia coli O157:H7. Epidemiological analysis of enhanced surveillance questionnaire data identified handling raw beef and shopping from the same national retailer (retailer A) as the common exposure. Concurrently, a microbiological survey of minced beef at retail identified the same strain in a sample of minced beef sold by retailer A, providing microbiological evidence of the link. Between September and November 2019, a further four primary and two secondary cases infected with the same strain were identified; two cases developed haemolytic uraemic syndrome. None of the four primary cases reported consumption of beef from retailer A and the transmission route of these subsequent cases was not identified, although all four primary cases visited the same petting farm. Generally, outbreaks of STEC O157:H7 in the UK appear to be distinct, short-lived events; however, on-going transmission linked to contaminated food, animals or environmental exposures and person-to-person contact do occur. Although outbreaks of STEC caused by contaminated fresh produce are increasingly common, undercooked meat products remain a risk of infection.

The spatio-temporal distribution of COVID-19 infection in England between January and June 2020.

The spatio-temporal dynamics of an outbreak provide important insights to help direct public health resources intended to control transmission. They also provide a focus for detailed epidemiological studies and allow the timing and impact of interventions to be assessed.A common approach is to aggregate case data to administrative regions. Whilst providing a good visual impression of change over space, this method masks spatial variation and assumes that disease risk is constant across space. Risk factors for COVID-19 (e.g. population density, deprivation and ethnicity) vary from place to place across England so it follows that risk will also vary spatially. Kernel density estimation compares the spatial distribution of cases relative to the underlying population, unfettered by arbitrary geographical boundaries, to produce a continuous estimate of spatially varying risk.Using test results from healthcare settings in England (Pillar 1 of the UK Government testing strategy) and freely available methods and software, we estimated the spatial and spatio-temporal risk of COVID-19 infection across England for the first 6 months of 2020. Widespread transmission was underway when partial lockdown measures were introduced on 23 March 2020 and the greatest risk erred towards large urban areas. The rapid growth phase of the outbreak coincided with multiple introductions to England from the European mainland. The spatio-temporal risk was highly labile throughout.In terms of controlling transmission, the most important practical application of our results is the accurate identification of areas within regions that may require tailored intervention strategies. We recommend that this approach is absorbed into routine surveillance outputs in England. Further risk characterisation using widespread community testing (Pillar 2) data is needed as is the increased use of predictive spatial models at fine spatial scales.

Utility of whole-genome sequencing during an investigation of multiple foodborne outbreaks of Shigella sonnei.

In April 2018, Public Health England was notified of cases of Shigella sonnei who had eaten food from three different catering outlets in England. The outbreaks were initially investigated as separate events, but whole-genome sequencing (WGS) showed they were caused by the same strain. The investigation included analyses of epidemiological data, the food chain and microbiological examination of food samples. WGS was used to determine the phylogenetic relatedness and antimicrobial resistance profile of the outbreak strain. Ultimately, 33 cases were linked to this outbreak; the majority had eaten food from seven outlets specialising in Indian or Middle Eastern cuisine. Five outlets were linked to two or more cases, all of which used fresh coriander although a shared supplier was not identified. An investigation at one of the venues recorded that 86% of cases reported eating dishes with coriander as an ingredient or garnish. Four cases were admitted to hospital and one had evidence of treatment failure with ciprofloxacin. Phylogenetic analysis showed that the outbreak strain was part of a wider multidrug-resistant clade associated with travel to Pakistan. Poor hygiene practices during cultivation, distribution or preparation of fresh produce are likely contributing factors.

Application of kernel smoothing to estimate the spatio-temporal variation in risk of STEC O157 in England.

Identifying geographical areas with significantly higher or lower rates of infectious diseases can provide important aetiological clues to inform the development of public health policy and interventions designed to reduce morbidity. We applied kernel smoothing to estimate the spatial and spatio-temporal variation in risk of STEC O157 infection in England between 2009 and 2015, and to explore differences between the residential locations of cases reporting travel and those not reporting travel. We provide evidence that the distribution of STEC O157 infection in England is non-uniform with respect to the distribution of the at-risk population; that the spatial distribution of the three main genetic lineages infecting humans (I, II and I/II) differs significantly and that the spatio-temporal risk is highly dynamic. Our results also indicate that cases of STEC O157 reporting travel within or outside the UK are more likely to live in the south/south-east of the country, meaning that their residential location may not reflect the location of exposure that led to their infection. We suggest that the observed variation in risk reflects exposure to sources of STEC O157 that are geographically prescribed. These differences may be related to a combination of changes in the strains circulating in the ruminant reservoir, animal movements (livestock, birds or wildlife) or the behavior of individuals prior to infection. Further work to identify the importance of behaviours and exposures reported by cases relative to residential location is needed.

Epidemiology of Hepatitis E in England and Wales: A 10-Year Retrospective Surveillance Study, 2008-2017.

Indigenous, foodborne transmission of hepatitis E virus genotype 3 (HEV G3) has become recognized as an emerging problem in industrialized countries. Although mostly asymptomatic, HEV G3 infection has a range of outcomes, including mild illness, severe acute hepatitis, and, of particular concern, chronic progressive hepatitis in immunocompromised patients. Public Health England has monitored cases of acute HEV infection in England and Wales since 2003. Between 2010 and 2017, enhanced surveillance using 2 linked laboratory databases and questionnaires on clinical features and risk factors was conducted. There was a year-on-year increase in the number of infections from 2008 (183) through 2016 (1243). Then, in 2017, the number of infections declined (to 912). As reported previously, HEV G3 group 2 (also known as "G3 abcdhij") is the predominant cause of acute infections, and older men are most at risk. Consumption of pork and pork products was significantly higher among patients than in the general population, but other previously reported associations, such as consumption of shellfish, were not observed. Ongoing surveillance is required to monitor future trends and changes in the epidemiology of the virus. The changing methods of animal husbandry and processing and distribution of animal products needs to be further investigated.

Analysis of the Cryptosporidium spp. and gp60 subtypes linked to human outbreaks of cryptosporidiosis in England and Wales, 2009 to 2017.

BACKGROUND: Cryptosporidium spp. are important causes of gastroenteritis that can be transmitted from humans and animals. We elucidated the distribution of species and gp60 subtypes in human outbreaks classified by transmission vehicle. METHODS: We used a combined database of national outbreak surveillance and reference unit data to analyse outbreaks by setting, vehicle, season, and linkage with suspected sources. RESULTS: A total of 178 outbreaks involving 4031 laboratory confirmed cases were identified; 82 (46%) outbreaks involved recreational waters, 74 (42%) animal contact, 4 (2%) environmental contact, 4 (2%) person-to-person spread, 3 (2%) food, 2 (1%) drinking water supplies, and 9 (5%) were of unknown source. The infecting Cryptosporidium sp. was identified in 131 (74%) outbreaks; 69 were C. parvum, 60 C. hominis, and in two outbreaks cases were infected with either species. Animal contact, environmental contact, and food-borne outbreaks were exclusively C. parvum and were mainly in first half of the year. Recreational water outbreaks were predominantly C. hominis and were mainly in the second half of the year. Outbreaks attributed to person-to-person spread were exclusively C. hominis and all occurred in October. Both C. parvum and C. hominis caused drinking waterborne outbreaks. Gp60 subtypes were identified from patients in 48 C. parvum and 38 C. hominis outbreaks, revealing more subtypes among C. parvum (n = 14) than C. hominis (n = 7) outbreaks. Cryptosporidium hominis IbA10G2 predominated (30 outbreaks). Of C. parvum subtypes, IIaA15G2R1 predominated (17 outbreaks), followed by IIaA17G1R1 (12 outbreaks), IIaA19G1R1 (four outbreaks), and other subtypes caused three or fewer outbreaks each. Linkage between cases and suspected sources by gp60 subtype was established in nine animal contact, three swimming pool, and one drinking water outbreak. CONCLUSIONS: The public health benefit of identifying infecting species and subtypes was twofold: (i) identifying and strengthening epidemiologic links between cases; and (ii) indicating possible exposures and sources to inform outbreak management. Gp60 subtype refined the epidemiological investigations, but a multilocus genotyping scheme would provide further benefit. Characterisation of Cryptosporidium spp. and subtypes needs to shift from predominantly supporting outbreak investigations to becoming nationally systematic.

Seasonality and the effects of weather on Campylobacter infections.

BACKGROUND: Campylobacteriosis is a major public health concern. The weather factors that influence spatial and seasonal distributions are not fully understood. METHODS: To investigate the impacts of temperature and rainfall on Campylobacter infections in England and Wales, cases of Campylobacter were linked to local temperature and rainfall at laboratory postcodes in the 30 days before the specimen date. Methods for investigation included a comparative conditional incidence, wavelet, clustering, and time series analyses. RESULTS: The increase of Campylobacter infections in the late spring was significantly linked to temperature two weeks before, with an increase in conditional incidence of 0.175 cases per 100,000 per week for weeks 17 to 24; the relationship to temperature was not linear. Generalized structural time series model revealed that changes in temperature accounted for 33.3% of the expected cases of Campylobacteriosis, with an indication of the direction and relevant temperature range. Wavelet analysis showed a strong annual cycle with additional harmonics at four and six months. Cluster analysis showed three clusters of seasonality with geographic similarities representing metropolitan, rural, and other areas. CONCLUSIONS: The association of Campylobacteriosis with temperature is likely to be indirect. High-resolution spatial temporal linkage of weather parameters and cases is important in improving weather associations with infectious diseases. The primary driver of Campylobacter incidence remains to be determined; other avenues, such as insect contamination of chicken flocks through poor biosecurity should be explored.

Utility of Whole Genome Sequencing To Describe the Persistence and Evolution of Listeria monocytogenes Strains within Crabmeat Processing Environments Linked to Two Outbreaks of Listeriosis.

This article describes the identification and investigation of two extended outbreaks of listeriosis in which crabmeat was identified as the vehicle of infection. Comparing contemporary and retrospective typing data of Listeria monocytogenes isolates from clinical cases and from food and food processing environments allowed the detection of cases going back several years. This information, combined with the analysis of routinely collected enhanced surveillance data, helped to direct the investigation and identify the vehicle of infection. Retrospective whole genome sequencing (WGS) analysis of isolates provided robust microbiological evidence of links between cases, foods, and the environments in which they were produced and demonstrated that for some cases and foods, identified by fluorescent amplified fragment length polymorphism, the molecular typing method in routine use at the time, were not part of the outbreak. WGS analysis also showed that the strains causing illness had persisted in two food production environments for many years and in one producer had evolved into two strains over a period of around 8 years. This article demonstrates the value of reviewing L. monocytogenes typing data from clinical cases together with that from foods as a means of identifying potential vehicles and sources of infection in outbreaks of listeriosis. It illustrates the importance of reviewing retrospective L. monocytogenes typing alongside enhanced surveillance data to characterize extended outbreaks and inform control measures. Also, this article highlights the advantages of WGS analysis for strain discrimination and clarification of evolutionary relationships that refine outbreak investigations and improve our understanding of L. monocytogenes in the food chain.

A spatial and temporal analysis of risk factors associated with sporadic Shiga toxin-producing Escherichia coli O157 infection in England between 2009 and 2015.

Infection with STEC O157 is relatively rare but has potentially serious sequelae, particularly for children. Large outbreaks have prompted considerable efforts designed to reduce transmission primarily from food and direct animal contact. Despite these interventions, numbers of infections have remained constant for many years and the mechanisms leading to many sporadic infections remain unclear.Here, we show that two-thirds of all cases reported in England between 2009 and 2015 were sporadic. Crude rates of infection differed geographically and were highest in rural areas during the summer months. Living in rural areas with high densities of cattle, sheep or pigs and those served by private water supplies were associated with increased risk. Living in an area of lower deprivation contributed to increased risk but this appeared to be associated with reported travel abroad. Fresh water coverage and residential proximity to the coast were not risk factors.To reduce the overall burden of infection in England, interventions designed to reduce the number of sporadic infections with STEC should focus on the residents of rural areas with high densities of livestock and the effective management of non-municipal water supplies. The role of sheep as a reservoir and potential source of infection in humans should not be overlooked.

National outbreak of Shiga toxin-producing Escherichia coli O157:H7 linked to mixed salad leaves, United Kingdom, 2016.

We investigated a large outbreak of Escherichia coli O157 in the United Kingdom (UK) with 165 cases between 31 May and 29 July 2016. No linked cases were reported in other countries. Cases were predominately female (n = 128) and adult (n = 150), 66 attended hospital and nine had features of haemorrhagic uraemic syndrome. A series of epidemiological studies (case-control, case-case, ingredients-based and venue-based studies) and supply chain investigations implicated mixed salad leaves from Supplier A as the likely outbreak vehicle. Whole genome sequencing (WGS) indicated a link with strains from the Mediterranean and informed the outbreak control team to request that Supplier A cease distributing salad leaves imported from Italy. Microbiological tests of samples of salad leaves from Supplier A were negative. We were unable to confirm the source of contamination or the contaminated constituent leaf although our evidence pointed to red batavia received from Italy as the most likely vehicle. Variations in Shiga toxin-producing E.coli surveillance and diagnosis may have prevented detection of cases outside the UK and highlights a need for greater standardisation. WGS was useful in targeting investigations, but greater coverage across Europe is needed to maximise its potential.

Whole genome sequencing reveals an outbreak of Salmonella Enteritidis associated with reptile feeder mice in the United Kingdom, 2012-2015.

Analysis of whole genome sequencing data uncovered a previously undetected outbreak of Salmonella Enteritidis that had been on-going for four years. Cases were resident in all countries of the United Kingdom and 40% of the cases were aged less than 11 years old. Initial investigations revealed that 30% of cases reported exposure to pet snakes. A case-control study was designed to test the hypothesis that exposure to reptiles or their feed were risk factors. A robust case-definition, based on the single nucleotide polymorphism (SNP) profile, increased the power of the analytical study. Following univariable and multivariable analysis, exposure to snakes was the only variable independently associated with infection (Odds ratio 810 95% CI (85-7715) p < 0.001). Isolates of S. Enteritidis belonging to the outbreak profile were recovered from reptile feeder mice sampled at the retail and wholesale level. Control measures included improved public health messaging at point of sale, press releases and engagement with public health and veterinary counterparts across Europe. Mice destined to be fed to reptiles are not regarded as pet food and are not routinely tested for pathogenic bacteria. Routine microbiological testing to ensure feeder mice are free from Salmonella is recommended.

Cyclosporiasis in travellers returning to the United Kingdom from Mexico in 2017: lessons from the recent past to inform the future.

During the summers of 2015 and 2016, the United Kingdom experienced large outbreaks of cyclosporiasis in travellers returning from Mexico. As the source of the outbreaks was not identified, there is the potential for a similar outbreak to occur in 2017; indeed 78 cases had already been reported as at 27 July 2017. Early communication and international collaboration is essential to provide a better understanding of the source and extent of this recurring situation.

Challenges in developing methods for quantifying the effects of weather and climate on water-associated diseases: A systematic review.

Infectious diseases attributable to unsafe water supply, sanitation and hygiene (e.g. Cholera, Leptospirosis, Giardiasis) remain an important cause of morbidity and mortality, especially in low-income countries. Climate and weather factors are known to affect the transmission and distribution of infectious diseases and statistical and mathematical modelling are continuously developing to investigate the impact of weather and climate on water-associated diseases. There have been little critical analyses of the methodological approaches. Our objective is to review and summarize statistical and modelling methods used to investigate the effects of weather and climate on infectious diseases associated with water, in order to identify limitations and knowledge gaps in developing of new methods. We conducted a systematic review of English-language papers published from 2000 to 2015. Search terms included concepts related to water-associated diseases, weather and climate, statistical, epidemiological and modelling methods. We found 102 full text papers that met our criteria and were included in the analysis. The most commonly used methods were grouped in two clusters: process-based models (PBM) and time series and spatial epidemiology (TS-SE). In general, PBM methods were employed when the bio-physical mechanism of the pathogen under study was relatively well known (e.g. Vibrio cholerae); TS-SE tended to be used when the specific environmental mechanisms were unclear (e.g. Campylobacter). Important data and methodological challenges emerged, with implications for surveillance and control of water-associated infections. The most common limitations comprised: non-inclusion of key factors (e.g. biological mechanism, demographic heterogeneity, human behavior), reporting bias, poor data quality, and collinearity in exposures. Furthermore, the methods often did not distinguish among the multiple sources of time-lags (e.g. patient physiology, reporting bias, healthcare access) between environmental drivers/exposures and disease detection. Key areas of future research include: disentangling the complex effects of weather/climate on each exposure-health outcome pathway (e.g. person-to-person vs environment-to-person), and linking weather data to individual cases longitudinally.